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Abstract 

PsJ ' We study the asymptotics of solutions of the Boltzmann equation describing the 

kinetic limit of a lattice of classical interacting anharmonic oscillators. We prove 
^ ■ that, if the initial condition is a small perturbation of an equilibrium state, and 

. vanishes at infinity, the dynamics tends diffusively to equilibrium. The solution is 

Q , the sum of a local equilibrium state, associated to conserved quantities that diffuse 

I to zero, and fast variables that are slaved to the slow ones. This slaving implies the 

I Fourier law, which relates the induced currents to the gradients of the conserved 



quantities. 



o 

B 

> 

H ; 1 Introduction 

If a piece of solid is heated locally and left to cool, the initial temperature distribution 
will diffusively relax to a constant temperature. The process is accompanied by a heat 
flow that is proportional to the local temperature gradient, according to the Fourier law. 
A mathematical understanding of this phenomenon starting from a microscopic model of 
matter is a considerable challenge [2]. 

A simple classical mechanical system modeling heat transport in solids is given by a 
system of coupled oscillators organized on a dimensional cubic lattice Z"^. The oscillators 
are indexed by lattice points x G Z*^, and carry momenta and coordinates {px, Qx)- In the 
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simplest model px and qx take real values and their dynamics is generated by a Hamiltonian 
which is a perturbation of a harmonic system: 

H{q,p) = \Y.p'- + \Y. 1-1y^''(^ - ^) + ^ I] (1-1) 

X xy X 

In ( 11.11) , the oscillators are coupled by harmonic forces generated by the coupling matrix u"^ 
which is taken short range. The parameter A describing the strength of the anharmonicity 
is assumed to be small. The classical dynamics is given by Hamilton's equations 

dH 

Qx = Px, Px = -T^- (1-2) 
OQx 

The dynamics fll.2l) preserves the Gibbs measures in the phase space, formally given 

by 

Z-^e-'^^^'^'P^dqdp (1.3) 

and describing an equilibrium state with constant inverse temperature j3 = 1/T. One 
expects initial states that agree with (11. 3p "at infinity" to be attracted by (II. 3p under the 

flow ira . 

While such a result is well beyond current techniques, one approach is to study this 
problem in an appropriate small coupling limit, the kinetic limit. In that limit, one takes 

= Re, one rescales space and time by e and, in the limit as e — > 0, one formally arrives 
at an evolution equation for the covariance of a Gaussian measure of the form [8]: 

W{x, k, t) + —Vuo{k) ■ VxW{x, k, t) = RC{W){x, k, t), (1.4) 

where x G M'^, k eT'^ = jlixl^ . The relation of W to the microscopic model is 

/ e~''^yW[x, k,t)dk = \im(^Uj{y - z){qx/e+zqx/e-z)e + i{qx/e+yPx/e-y)e), (1-5) 

where (— )e is the state at time t/e. W is real and positive. Since the function W{Rx, k, Rt) 
satisfies (11.41) with i? = 1 we may, with no loss of generality, make that assumption. 
Equation (ll.4p is a Boltzmann type equation where the collision term is given by 



C{W){xXt) = ^ 



/ dkidk2dkz{(jJuJiUJ20J'i) ^5{uj + uji — uj2 — uJ^) 
6{k + ki-k2- k3)[WiW2W3 - W{WiW2 + W1W3 - W2W3)] (1.6) 



with uj = uj{k), W = W{x,k,t), uJi = Lj{ki), Wi = W{x,ki,t), i = 1,2,3. The sum 
fc + fci — ^2 — ^3 is mod 27rZ. We normalize J^a dk = 1. 
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Analogously to the collision term of the standard Boltzmann equation for gases, the 
collision term (11.61) describes a scattering process. The scattered "particles" are phonons, 
i.e. vibrational modes, and the scattering process in (11.61) involves four phonons. The 
phonons carry momenta, k,ki,k2,k3 and energies u{k),uj{ki),uj{k2),uj{k3) and the scat- 
tering process conserves the total energy and momentum. 

We may pose the problem stated in the beginning, i.e. the question of relaxation to 
equilibrium, for the equation (II. 4p . Thus, starting from some initial W that coincides "at 
infinity" with an equilibrium state we would like to prove that W{t) tends diffusively to 
equilibrium and the process is accompanied by heat fluxes governed by Fourier's law. In 
this paper, we prove such a theorem (Theorem 1 in Section 3), for an initial condition that 
is close to equilibrium. The equation (ll.4p becomes a coupled system for "fast" and "slow" 
variables (see (13. 6p . (13. 7p below). The slow variables correspond to the temperature and 
the chemical potential of another conserved quantity, defined in (I3.16p . The slow variables 
flow under a nonlinear diffusion equation whereas the fast ones, that include the heat 
currents, are " slaved" to the slow ones by the Fourier law, which is discussed after stating 
Theorem 1 in Section 3. 

There is a rich literature on the derivation of hydrodynamics from the Boltzmann 
equation for gases. In particular in [I], it is proven that, in an appropriate limit, the 
Boltzmann equation gives rise to the Navier-Stokes equations (see also [HI El [10] for more 
results and background on the link between Boltzmann and Navier-Stokes). In Theorem 
2 (Section 3), we consider the diffusive scaling limit, i.e. the limit where one scales the 
spatial variable by e, and the time variable by e^, and we show that a suitably rescaled 
solution of the Boltzmann equation (11.41) solves, in the limit e — > 0, a nonlinear heat 
equation. 

Compared with the Boltzmann equation for gases, eq. (11.41) has some nice features. 
In particular, its long time dynamics is described by hydrodynamic equations that are 
considerably simpler than the Navier-Stokes- Fourier system for the former (in the Euler 
scahng, where space and time are both scaled by e, eq. (I1.4p reduces to the Fourier 
law). A mathematical theory of (ll.4p however appears to be lacking. This paper should 
be considered as a step in that direction. For an attempt to prove the Fourier law for 
phonons, beyond the Boltzmann approximation, see [21 H]. 

2 The collision term 

The collision term (11.60 involves integration over the subset determined by the constraints 
in the delta functions. Some assumptions on u are needed for regularity of the resulting 
measure. The standard nearest neighbor coupling between the oscillators corresponds to 
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UJ = uj^ with 

d 

a;°(A;)2 = 2^(1- COS A;^) + r (2.1) 
i=i 

where r > is the pinning parameter. The regularity properties stated in Proposition 2.1 
should hold for this u as long as r > 0. The proof of this seems, however, quite tedious 
and, therefore, we take 

uj{k) = uj'^{k)\ (2.2) 

with r > 0, which is simpler to analyze. Moreover, we need to take the spatial dimension 
d>2. From now on, we'll make that assumption. 

For r = and u!{k) = uj^{k) the regularity problems are more serious and we expect 
to need d > 3. Moreover, the spectrum of the linearized collision operator has no gap 
and the nature of the flow of eq. ( 11. 4p is quite different since also the slaved modes are 
diffusive. Progress in that case is so far hampered by the complexity of the integrals 
appearing in (11. 6p . 

The basic properties of the collision term are summarized in Proposition 2.1. Since 
the X variable plays no role in C we consider C as a map in 

B = C°(T'^) 

equipped with the sup norm, denoted by || • ||. We also need the Hilbert space 

n = L\T^,uj\k)dk), (2.3) 

whose norm will be denoted by || ■ Ht^. In this paper, we shall use C or c to denote 
constants that may vary from place to place. 

Proposition 2.1 (a) C : B ^ B with \\C{W)\\ < C\\Wf, for some constant C < oo. 

(b) The equation C(W) = has, for W >0, exactly a two parameter family of solutions 
in B 

where W > for A > — r^. 

(c) For allW eB 

' dkC{W)ik) = J dkuj{k)C{W)ik) = 0. (2.5) 

(d) Let —L be the linearization of C{W) around Wq = uj{k)~^ i.e. 

L = -DC{Wo). (2.6) 
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Then, L is bounded in B and hounded and positive inTi. It has two zero modes 

Luj-^ = Lcu"^ = 0. (2.7) 

and the rest of its spectrum in B and Ti is contained in 3fJA > a > 0. 

Moreover L = M + K where M is a multiplication operator by a strictly positive 
continuous function and K is an integral operator compact in Ti and B and satisfying 



sup / \K{k,k')\dk' < oo. (2i 

k J 



Proof, (b) We use an argument given in [8], for a similar kernel. First, writing the [— ] 
in (11. 6p as 

3 

J] W{h) [W{ko)-' - W{ks)-' - W{k2)-' + W{k,r'] , (2.9) 

i=0 

we have: 

3 

CiW)ik) = ^ / f\idhu~'Wih))5iuJo + CU,-U2- CUs) 

■6iko + k,-k2- k,) [Wiko)-' - Wik,)-' - Wik^y + Wik,)-'] , (2.10) 

which vanishes, for W{ki) > 0, if f{k) = W{k)~^ is a collisional invariant, i.e. a function 
satisfying 

f{ko) + f{k,) = f{k2) + fiks), (2.11) 

on the set of vectors ki in the support of the delta functions in (12.101) . In [9] it is proven 
that all L^-solutions of (12.111) are of the form au{k) + b. Hence, the functions in (12.41) 
satisfy C{W) = 0. 

To show that these functions are the only solutions of C{W) = 0, write, for g bounded. 



[ dkg{k)C{W){k) = ^ [ Y\{dkiU-'Wih))5{iUo + 

J 4 Jjid fj- 



UJl — UJ2 — UJ^) 

i=0 

■6{ko + ki-k2- ks)g{ko) [W{ko)-' - W{ks)-' - W{k2)-' + W{ki)-'] , (2.12) 

and use the symmetry ^ 1 of the rest of the integrand, to replace g{ko) in front of 
[— ] in (12.121) by ^{g{ko) + g{ki), and, using the antisymmetry of the rest of the integrand 
under the exchange ^ 3, 1 ^ 2, we may replace h{g{ko) + g{ki) by — i(fi'(^3) + fl'l^a))- 
This implies that (I2.12p is proportional to 

/ TT('^^it^i ^W{ki))6{uo + UJ1-UJ2- uj3)6{ko + ki- k2- k^) 
J-^'" f=o 

■ [g{k,) + g{k,) - g{k2) - g{k,)] [W{k,)-^ - W^k,)'^ - W{k2)-^ + W{k{)-'] (2.13) 
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Taking g{k) = W{k) ^, we obtain 
h^^ f=o 

■ [w{ko)-' - wik^r' - wik^r' + w{k,r'] ' , (2.14) 

which vanishes, for W{ki) > 0, only if 14^ (A;) is a colhsional invariant. 

(c) Since / dkg{k)C{W){k) is proportional to f l2.13p . it vanishes, for any W, whenever g 
is a collisional invariant, in particular for g{k) = 1 and g{k) = uj{k). 

(d) Differentiating C{Wt,a) = with respect to T and A at T = 1 and A = 0, we obtain 

Luj-^ = Luj-^ = 0. (2.15) 
Explicitely, we have (note the — sign in (12. 6p ) 

Lf = / TT dkiUJ~'^8{ujQ + UJ1-UJ2- i^3)S{ko + ki-k2- k^) 

{uikoffiko) + uj{k,ff{k,) - U\k2)f{k2) -uj\k,)f{k^)) . (2.16) 

To obtain (I2.16P from (11.61) . write the [— ] in (11.61) as in (12.91) and expand 

W{kir' = {u{h)-^ + f{h)Y^ = u{h) - u{hff{h) + ... (2.17) 

If we take the first term in (12.17^ in the [— ], we get because of the delta function in (II. 6p . 
So, the only linear terms in / correspond to taking the second term in (I2.17P inside the 
[— ], and replacing IlLo ^(^«) IlLo'^j"^ which implies (I2.16p . Using the symmetries 
that led to (I2.14p . we see that (/, Lf) is proportional to 

/ TT dkiUi'^5{iji!o + uji- UJ2- oj-i)5{kQ + ki - k2 - k^) 

{uikoffiko) + Uj{hff{h) - Uj\k2)f{k2) -uj\k^)f{k,)f > 0. (2.18) 



Thus, the zero modes of L are of the form / = uj^^g, where g are collisional invariants, 
i.e. of the form auo{k) + h (see part (b)) and therefore (12.71) give the only zero modes in 
H, and in B. 

From (12.160 , we get the decomposition L = M + K with 
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4^0 



/ ITrf/cjCUj ^(5(u;o + cui - u;2 - u;3)(5(/co + A^i - /i;2 - ^3)t^^(^o), (2-19) 



where k = kQ, and K is an integral operator 



{Kf){k)= / K{k,k')f{k')u{k'fdk' 



given by the last three terms in fl2.16p . i.e. 



K{k, k') = -^{uj{k)uj{k'))-\l'^^\k, k') + A;')) (2.20) 



where 

I^^\k,k') 



2 / dki{uj{ki)uj{ki + k-k')Y^5{uj{ki)-uj{ki + k-k')+uj{k)-uj{k')), (2.21) 

and 

I^^\k,k') = 

dki{uj{ki)uj{ki + k- k')y^6{-uj{ki) - uo{ki + k-k')+ uj{k) + uj{k')). {2. 22) 

The singularities of J^^-* and I^"^^ are similar, and we will study J*-^-* first; let 

r='-{k-k') 

and change variables in (12.211) to ki = g — r — (|, . . . , |) and write | for (|, . . . , |). 
Recalling ^ we get 

I^'\k,k')= [ dq6iniq,r,k')) giq,r) (2.23) 

with g = {uj{q — r — ^■n))uj{q + r — ^tt)))^^, and 

d 

n{q,r, k') = ^(sin(gj + rj) - sm{qj - rj)) + uj{k' + 2r) - u{k'). (2.24) 
i=i 

The summand equals 2 cos sin . We change variables to cos qj = Sj{l — yj) with Sj = 1 
in the region qj G [— |, |], and Sj = —1 in the region \qj\ > |. In both cases, yj G [0, 1]. 
Our integral becomes a sum of integrals of the form 

Is{k,k')= [ 6{2y2yjSjsmrj-m{k',r))G{y,r)l\y;h{y,)dyj, (2.25) 

with h{y) = (2 — y)~^ , G smooth, and 

m{k\r) =uj{k' + 2r) - u{k') +2^SjSmrj. (2.26) 

3 

Let sinrj ^ for all j i.e. /cj — k'^^Q modulo 27r and m(A;', r) 7^ 0. Given fc, the set of k' 
such that these conditions hold is of full measure. Then, 

Is{k,k') = \[\smr,\-'^ Js{k,k'), (2.27) 
3 
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where 

Js{Kk') = / J](^/,(^^))5(2^y,.s,-m(A:',r))G(y„r), (2.28) 

with = yj/\ sinrjl and = ±1. 

Consider e.g. d = 2. Integrating the delta function over the variable y2 shows that 
the singularities of J are the same as those of the function 

/"-'x(^. + r»)e|o.|smr.||)% 

Jo ^/yVsy + rn 

evaluated at m = s'm{k', r) where s, s' = ±1. j is bounded by 

\j{r,k',m) \ < C(l + |log|m||). (2.30) 

Since m(A;', r) is real analytic, Y[j | sin r^p^ log m is square integrable and thus / G L^(T'^x 
T'^) and, so, K, given by fl2.20p . is compact in H. From fl2.16p . we see that L is self adjoint 
in n. 

To show compactness in B, we will use Holder continuity properties of Is{k,k'). We 
want to bound 

Jh {k,k')-h{k,k')\ 



sup / dk 



for some a > 0. We will analyze the dependence in k of each factor in fl2.27p . 02.281) 
separately. For m, f G (0, 1] we have, for any n, 



I -1 -1| I — — I I " + l-k J_ n + 1 n + 1 

\u ^ — V ^ I = |-U2" — tiZ" 1 1 u '^"V 2" I < <-^n|w — f 1 2" maxj-u 2™ , f 2™ | (2.31) 



k=l 



smce U2n can use this, with n > 1, to control the k dependence 

of the factor Ylj |sinrj|~^ in fl2.27p . For Jg, let us, for simplicity, again study j. Let 
rh = m{k',f), with f = |(A; — k'). Applying (12.311) to (I2.29p we get, with o: = 

/•|sinri| -j^ ^ 

< Ca|m — m|"(l + max{|m|^", |m|^"}) 

< Ca\k - kl^il + max{|m|-°, \m\~"}), (2.32) 

where the last inequality holds because m is analytic in k. Next, let | sinfi| = | sinri| +6, 
assuming S > 0. Then, using the fact that 

f\^^-r,\+^ ^(sy + m)e[0,\^nr,\])dy , S 

/ r- I I ,< Cmm(l + |log|m||,log(l + — r)), 

ilsinr-il ^^sy + m |smri| 



we get the following crude estimate: 

1^2, k' , m) — j {ri, r2, k' , m)\ < C(52 + x(| sinri| < (1 + | log |m| |)). (2.33) 

A similar estimate holds for the dependence with respect to r2, and for the one with respect 
tommxisy + m) G [0, | sinr2|]). Since m is real analytic, Hj | sin rj|~^|m(/c', ^ (A; — fc'))|~" 
is integrable in k' uniformly in k for a > small enough. Moreover, the integral of 
Y[j I sin rj\~^\m{k', ^{k — k'))\~°'x{\ sinri| < 62) over k[ is of order 5", for a small. 
We get, by combining all these estimates, that for some a > and all f E B, 



Is{k,k')f{k')dk'\\^<C\\f\\ (2.34) 

where || — is the C° norm. We may repeat this analysis for writing: 

I^^\k,k' + n) = 



- I dki{uj{ki - n)iu{ki + k- k')y^5{uj{ki) - uj{ki + k - k') + u{k) - cj(fc')I^.35) 



where we changed variables ki ki + [it , . . . , tt) in the integrand and used formula (12. ip . 
(12. 2p for uj{k) and cos(A;j + tt) = — cos kj. Since the factor {uj{ki — 7i)u{ki + k — k'))~'^ is 
smooth, we can do exactly the same analysis as above and obtain the bound (12.340 for 
Then, the same bound holds for the kernel given by (l2.20p .This shows that K is 
a continuous map from B into the space of Holder continuous functions. Since C" is 
compactly embedded in B (by Ascoli's theorem), this proves that K is compact. Since, 
using (E2Q]), (E2ID, (E22]), (EJSI), 



M{k) = / dk'K{k,k')uj{k 



'\2 



we also get M & B (actually M is Holder continuous) and that L is bounded (since K is 
compact). 

The case ci > 2 is similar, 
(a) We proceed as in (d) with the delta functions and get 

\CiW)ik)\<C\\W\\' [ dk' [ 6i2Ty^s,sinr^-mik',r))l\y;h{y^)dy, (2.36) 

which, using (I2.27p . (12.290 . proves the claim. □ 
The identities (12.51) yield two conservation laws. To show this, define, for a = 1,2, 

Ux,t) = -^{u;-'^,VuW), (2.37) 
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where the scalar product is in Ti, see (12. 3p : ji is the thermal current and j2 can be called 
the phonon number current. Similarly, set 

T„(x,t) = (2.38) 

Ti is the temperature and T2 is related to the phonon chemical potential. Equations (11.41) 
and (12. 5p give then the conservation laws 

T; = V-j„. (2.39) 

Tq are the slow modes that will diffuse. The currents ja will be related to their gradients 
via the Fourier law, see (13.171) below. 



3 Results 

Let W{p,k,t) denote the Fourier transform of W in the x variable. We shall look for 
solutions of (11.41) of the form: 

W{p, k, t) = W,{k)5{p) + w{p, k, t). (3.1) 

where, we recall, Wq = u^^. The equation (II. 4p becomes then, 

w = —Dw + n{w), (3.2) 

where the linear operator is given by 

D = L + ^p-Vuj{k) (3.3) 

and the nonlinear term is n{w) = C{Wo + w) + Lw = C(Wo + w) — DCiWo)w. Written 
as an integral equation, (13. 2p becomes 

w{t) = e-*^u;(0) + r dse-^'-'^^'niwis)). (3.4) 







We need to decompose this in terms of the slow and the fast variables. Let P be the 
orthogonal projection in the Hilbert space Ti on E = span {uj^^^uj""^} and let Q be the 
one on the complement of E. The identities (12. 5p can then be written as: 

Pn = 0, (3.5) 

or n = Qn, since, by differentiation, (12. 5p implies the same identities with C{W) replaced 
by DC{Wo)w. Let 

Pw = T, Qw = V. 
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Then, w{t) = T(t) + v{t), and eq. (13.21) becomes 

T{t) = Pe-'^w{0) + / dsPe-'^'-'^^Qn{w{s)) (3.6) 

Jo 



v{t) = Qe-'^w{0) + / dsQe~^'-'^^Qn{wis)) (3.7) 

Jo 

Define 

e{p,t):={l + {t + l)p'r^ (3.8) 

with n > d/2 and let St be the space of continuous functions f{p, k) equipped with the 
norm 

ll/lli^sup eip,tr'\fip,k)\=snp e{p,t)-'\\f{p,-)\\ (3.9) 

k,p p 

and let S be the space of functions w{t) G £t, with 

ll^llf = sup ||w(t)||t. (3.10) 

t 

Let K, : E E he the linear operator 

K= {27v)-'^PdiUjL-^diUjP, (3.11) 

where di denotes the derivative with respect to the first argument, k, is strictly positive 
by Proposition 2.1.d, since the set diUjE forms a two dimensional subspace in E-^ (by 
symmetry: diuo is odd in A;, while the functions in E are even). Let 

ro(p,t) = e-*^''^T(p,0)x(N<l), (3.12) 

and 

vo{p. t) = ^L-'p ■ VuToip, t) (3.13) 

(these lie in E and E-^ respectively). Then, we have our main result: 

Theorem 1. There exists 6 > such that, for \\w{0)\\o < 6, the equation (13.21) has a 
unique solution in S, w{t) = T{t) + v{t), satisfying, for t>l, 

\\T{t)-T,{t)\\t<C(t)t-^\wmo, (3.14) 
\\v{t)-vomt<C{t)t-^\\w{mo: (3.15) 

where C{t) = Clog(l + t) for d = 2 and is constant for d > 2. 

Remark. The norm (13.91) and (I3.14p imply that, in x space, T{x,t) is given by the usual 
diffusive Gaussian term, which, for bounded x, decays like t~2 , plus a correction bounded. 
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uniformly in x, by 0(t-^). Eq. (I37[5|) says that the fast mode v is slaved to the slow one 
T, i.e. it decays like an explicit term, the Fourier transform of fo(p, t), plus a correction 
uniformly bounded by 0{C{t)t~^). 

Eq. f l3.15p also implies the Fourier law for the leading terms of the solution. Write T 
in the basis 

T{pXt)= Y,uj-I'{k)fp{p,t), (3.16) 

13=1,2 

(since the basis is not orthogonal, Tp does not coincide with the Fourier transform of Tp 
in fl2.38p ). Then, in x-space, the currents fl2.37p (where, by symmetry, only the v part of 
W contributes) become, up to terms 0{C{t)t ^), 

ja = /ta/jVT^, (3.17) 

P 

with the positive conductivity matrix 

Ko,p = {UJ''',KUJ-^), (3.18) 

since, by symmetry, ^diUjL~^djUJuj~^) equals zero for i ^ j, and, for i = j, i can be 
chosen equal to 1. 

Finally, we can also derive a nonlinear heat equation as the hydrodynamic scaling limit 
of the Boltzmann equation (11.40 . We scale 

W{x,k,t) = W{ex,k,eH), 

where W{x,k,t) satisfies ( 11.40 ). Thus, we obtain, for W, the equation: 

W{x, k, t) + {2TTe)-^Vuj{k) ■ VW{x, k, t) = e-'^C{W){x, k, t). (3.19) 

We shall solve it with initial data W^|t=o = uj{k)~^ + w{x, k, 0) and 

w{x, k, 0) = T(x, k, 0) + ev{x, k, 0), (3.20) 

with T(x, -,0) G E, vo{x,-,0) G E^, Vx. W{x,k,0) = W{ex, k,0) has spatial variations 
at scale e~^. 
We have then 

Theorem 2. There exists 5 > such that, for all e < 1 and ||T(0)||o, ||'y(0)||o < S, the 
equation i\3.19\) has a unique solution = + ev^ G S. Moreover, T^{t) T{t) and 
v'^{t) — > v{t) in St, for all t > 0, as e 0, where T and v are the unique solutions of 

DC{u-^ + T)v = {2'k)-^Vu} -VT (3.21) 
f = -{2tx)-^PVuj -Vv (3.22) 
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Remark. Since we have DC{ijJ ^ + T) instead of DC{uj ^) = —L in fl3.22l) . we get a 
nonlinear heat equation for T: 

t = V • (/C(r)VT) (3.23) 

with 

/C(T) = -{27T)-'^PdiUjDC{uj-^ + T)-^diUjP, (3.24) 

(remembering that L is positive, we see that DC is negative and that /C(T) is a positive 
matrix). One can show that the long-time behavior of fl3.23p is similar to the one in 
Theorem 1. 



4 Proofs 

We start by proving some results on the spectrum of the operator D defined in fl3.3p . 
Recall that, by Proposition 2.1 d, ^{a{L)\{0}) > a. 

Proposition 4.1 There exists po > 0, and b > 0, such that, both in B and in Ti, the 
following holds: 

a) // \p\ < Po then 

ct(D) = {Ai, As} U a 

with inf 3?(T > a/2 and 

K=p'fi. + 0{\pf) (4.1) 

where /Xj > are the eigenvalues of the operator k, defined in (13.111) . 

b) // IpI > Po then mi^a{D) > b. 

Proof, a) Let 

L{f3) = L + /3p 

where p = {2TT)^^\p\^^p-'Vu!{k) so that D = L{i\p\). p is a bounded operator in B and in H. 
Hence, for \P\ < po small enough, the spectrum of L{P) consists of two eigenvalues close 
to the origin, with the rest of the spectrum having real part larger than a/2. Moreover, 
L{(3) is self adjoint in Ti for (3 real and thus forms an analytic Kato family (see [7j). 
Therefore, there is a po such that, for |/?| < po, the two eigenvalues of C{P), Ai, A2, and 
the associated eigenfunctions ipi, are analytic in \P\ < po- To compute their Taylor series, 
use the decomposition w = T + v and project the equation (L + Pp)w = \w onto QH 
and PH. Then, using PpT = (since pT is odd in k), we get a pair of equations: 

Lv + (3Qp{T + v) = Xv, 
l3Ppv = AT, 
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whereby, since 3?((t(L)\{0}) > a, we get, for /? small and A = 

V = -(3{L - A + pQp)'^QpT = -pL-^QpT + 0{p^)T, 

where denotes a bound on the operator norm in B. Since QpT = pT and T = PT, 

we get 

-p^PpL-^pPT + 0{I3^)T = XT. 

The 2x2 matrix PpL^^pP is strictly positive, by Proposition 2.1 (d). The result follows 
by letting (3 = i\p\, and by observing that, since the functions in E are invariant under 
reflections and permutations of the coordinates of k, the matrix K,{p) = PpL^^pP equals 
fi;(ei) i.e. the one given in fl3.11l) . 

b) Write D = Ai + i\p\p + K a.s a sum of a multiplication and a compact operator. Since 
3?cr(A^ + z|p|p) > inf > we need only to show that all eigenvalues A satisfy ReA > b, 
if bl > Po- We decompose H = He ® Ho into even and odd subspaces. Let A be an 
eigenvalue. Then, since p is odd, and the operator L is parity-preserving, 

LWe + i\p\pWo = XWe 

Lwo + i\p\pwe = Xwo (4.2) 

Write X = X + iy and suppose that x < | (where a = inf cr(L \-Ho))- Then {L — A) {-h^ is 
invertible, and (14. 2p becomes 

iL + p^p{L-X)-'p)we = Xwe. (4.3) 

Let llifellH = 1- Then, taking the scalar product of (14. 3 p in H with We, and writing 
separately the real and imaginary parts, one gets: 

X = {we, {L + p^p{L - x){{L - xY + yyp)we) 

>^p'{we,p{{L-xr + y')-'pWe), (4.4) 

since L > and, as operators in Ho, {L — x){{L — x)^ + y"^)'^ > ^{{L — xY + y"^)^^- For 
the imaginary part, we get: 

y = yp^iwe, pi{L - xY + y^Y^pWe). (4.5) 

If y 7^ (14. 4p and (14. 5 p imply x > a/2, against our assumption. Thus ?/ = 0, i.e. all 
eigenvalues X oiV with ReA < | are real, and, by (14. 4p . positive. 

Let 

/(r,A)= inf [w , {L + r"^ p{L - X)-^ p)w) . 

\\w\\=l 

Since, for X < a/2, (L — A)^^ > c > 0, we get 

(w, (L + r'^p{L - XY^p)w) > {w, {L + cr'^p^)w). 
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To bound this from below, observe that both terms in L + cr^p^ are positive, that 
{w,Lw) > c||Qw|p, and that {Pw, p'^Pw) > c||Pw7|p, by inspection of the functions 
in PTi and of p. So, if > c'||w|p, we can use the first lower bound, while, if 

IIQwiP < c'llti'lp, i.e. Il-Ptyp > (1 — c')||w||^, we can use the second one, if c' << c, to 
get L + cr^p^ > c"r^ for r < Pq, and po small enough. So, l{p, \) > c"p^, for \p\ < Pq, 
and A < a/2. Besides, /(r. A) < l{r', A) if r < r'. Since by (14. 3 p A > l{\p\, A) we conclude 
A > c"pI, for \p\ > Pq. Thus, we can take b = min(|, c"pl) □ 

We are now ready to state the estimates for the heat kernels in equations (13. 6p and 
([32D- Let 

R{z) = (z-D)-'^ 

be the resolvent of the operator D and 7 be a circle of radius a/4 around the origin. Then 

P=<fR{z)^, (4.6) 

is a projection onto the 2-dimensional eigenspace of D introduced in Proposition 4.1 a. 
Let Q = 1 - P. We have: 

Proposition 4.2 For all t > we have 

||e-*^|| < C(e^^*P' + e-^*) (4.7) 

where \\ ■ \\ is the operator norm in B. Moreover, there exists po > so that for \p\ <po 

||Pe-*^Q|| + ||Qe-*^P|| < C\p\{e"''''' + e"^*) (4.8) 
ilQe-*^Q|| < C{p\-''^" + e-^*) (4.9) 
||e-*^Q|| < Ce-"' (4.10) 



Proof. Since D = L + ^p ■ Vu!{k), we have, by the resolvent expansion: 

P= Lz- L)-' ^+<({z- L)-'^p . Vu{z - L)-' ^ + 0{p') 



= P - —L-^p ■ VujP - —Pp ■ VujL-^Q + C(p2) 
27r 27r 

= P + AP + PBQ + 0{p^) (4.11) 

where, here and below, 0{p'^) denotes a bound on the norm of operators in and we 
define 

A = — -L-^p-Viu, B = — -p-VujL-^. (4.12) 

2tx 27r 
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Hence, for \p\ small, 

\\PQ\\ + \\QP\\ + \\PQ\\ + \\QP\\ <C\p\. (4.13) 

Now, write 

e"*^ = Pe-'^P + Qe-'^Q, (4.14) 

(since P, and Q project on invariant subspaces of D, we have e~^^P = Pe^'^^P, e^^^Q = 
Qe-'^Q), and 

e-'^Q= L-'^R{z) ^ (4.15) 

where the curve 7 goes around the part of the spectrum of D that lies in the complement 
of a ball of radius a/4, centered at the origin. Since L is bounded and \p\ small, the length 
of 7 is 0{1), and we get (14.101) . Then, the other claims follow, for \p\ small, from fl4.13p . 
( 14.141) and Proposition 4.1 a. 

To get (14. 7p for \p\ > pq we note that by Proposition 4.1, e~*-^ is given by the right 
hand side of (I4.15p . where the curve 7 can be chosen to be a rectangle in > b, with 
vertical sides of length const \p\, and horizontal ones of length 0(1) (since L is bounded). 
Hence, on the curve 7, we have |e~*^| < e~^*, but the length of the contour of integration 
is not bounded as |j9| 00. To control the integral, recall that D = M + i\p\p + K, and 
let 

Ro{z) = {z-M-t\p\p)-\ 
Then, by the resolvent formula, 

R = Rq + RqKRo + RqKRqKR. (4.16) 

The first term in (I4.16p . inserted in (14.150 . gives the multiplication operator e"**^*^'^*!^'^-', 
which satisfies the bound b) with c = inf M. For the second term, we use |e~*^| < e"*** in 
the integral and write 

I (f dk'dzRo{z,k)K{k,k')Ro{z,k')\ < [ dk'\K{k,k')\ (fdz- 1 (4.17) 

where zi = (M + i\p\p){k), Z2 = (M + i\p\p){k'). Since \z — Zi\ > c on 7, we have 

dz. 1 r < C, (4.18) 

\Z — Zi\\z — Z2\ 

uniformly in k,k'. Since, moreover, J \K{k, k')\dk' < C, uniformly in k (which holds by 
(12. Sp ). we get the bound in b) for the second term in (I4.16p . The third term in (14.160 is 
similar, since KR is bounded. □ 

Proof of Theorem 1. Let us write eq. (13.41) as 

w{t) = wg{t) + N{w,t) (4.19) 
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where Wi{t) is the solution of the hnear problem: 

We = e'^^w{0). (4.20) 

We will solve f l4.19p in the space £. 

The linear term is bounded, using (14.71) . and recalling the definition (13.81) of the weight 
function e, by 

WweWe < Csupe(p,t)-ie(p,0)(e-'=*P' + e-^*)||«;(0)||o < C||uKO)||o. (4.21) 

p,t 

Consider then the nonlinear term in eq. (13. 4p . By an easy extension of Proposition 
2.1 a (see ^M)), 

IHwismi < CiMMsrip) + ||^|||e(.)*=^(p)). (4.22) 
We need the simple estimate 

e{sr\p)<C{l + s)-h{p,s), (4.23) 

which follows from 

e(s)*^(p) < 2 / dq e{p - q,s)e{q,s) < Ce{p,s) / dqe{q,s), (4.24) 

J\p-Q\>i\p\ J 

(the first inequality holds because either |g| > ^\p\ or \p — q\ > h\v\)- The last integral 
yields the power of 1 + s since n > \d m. (13. Sp . The last term in ( 14.22^ is smaller, for 
\Me<l- 

Thus, the second term in equation (13.60 and the second term in (13.70 are bounded as 

by 

\\RN{w,t,p)\\<C\\w\\l f ds\\Re-^'-'^''Q\\{l + s)-ieis,p) (4.25) 

Jo 

where R equals P in (13. 6p and Q in (13. 7p . Proposition 4.2 imphes 

||i?e-(*-^^^g|| < C{{1 + t - s)-i"^e-i(*~^)''' + e'^^*"^)), (4.26) 

with m = l if i? = P and m = 2 if i? = Q (using |ppe-ip'(*-') <C{l+t- s)"*"", for 
t — s > 1). Bounding 

e-fp'(*-«)e(s,p) < Ce{t,p) 

and 

e-f(t-«)e(s,p) < Ce{t,p) 
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we may bound the remaining s-integral in (14.251) . for m = 1,2, i.e. for m < d, hj 
C(t)(l + t)-^™ where C{t) may be taken to be constant, except for d = 2, where C{t) = 
Clog(l + t). Hence we end up with 

\\RN{w,t)\\t < C{t){l + t)-i'^\\w\\l (4.27) 

In particular, we have 

\\Niw)\\e<C\\w\\l 

so that in (14.191) maps a ball of radius e in S into itself, for e small enough and is 
obviously a contraction there. The existence of a solution for eq. (14.191) then follows from 
(14.211) and the Banach fixed point theorem and, since m = 1 if i? = P, i.e. for (13. 6p . and 
m = 2 if R = Q i.e. for (13. 7p . we obtain the bounds 

\\Tit)-Temt<Cmi+t)-i\\wml (4.28) 

\\v{t)-vemt<cm+t)-'\\wmi (4.29) 

where we wrote 

Wi = Pwi + Qwi = Te + ve. (4.30) 

To conclude the proof of Theorem 1 we need to relate Tq, vq to Te, ve. For this we 
need to write the leading terms of Wi{t) more explicitely. Let us formulate this a little 
more generally, which will be useful also in the proof of Theorem 2. Let us denote 

K{t) = e-'P'^P (4.31) 

where k is given in eq, (13. lip . We have then the following Lemma, proven at the end: 

Lemma 4.3 Let \p\ < pq. Then the semigroup e~^^ can be written with respect to the 
decomposition E © E-^, as 



_tD _ I K{t) K{t)B \ ( 0{\p\e-^'P^ +p'^e-^') ^(p^e-^*?' + |p|e-^*^ 



^ ' AK{t) AK{t)B + R{t) J ^ yOip^e-^'P' + \p\e-^') Odppe-'^*^' 

where O is with respect to the operator norm in B, 

R{t) = QQe-'^^QQ (4.32) 

and the operators A and B are defined in equation (14.121) . 

Returning to the proof of Theorem 1 and recalling the definition (I3.12p of Tq and 
(14.201) . (I4.30p of Ti, we get from Lemma 4.3, for \p\ < Pq, 

||T,(t, p) - To(t, p) II < C\p\ {e^^P' + e-^')e{p, 0) ||^(0) ||o. (4.33) 
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Since vq given in eq. (13.131) equals, see (14.121) . vq = AK(t)T{0) for \p\ < po < I, we have, 
from Lemma 4.3, 

\\veit,p) - Vo{t,p)\\ < Cip'e-^^' + e-^*)e(p, 0)||w(0)||o (4.34) 
where we used the bound, coming from (14.101) . 

\\R{t)\\<Ce-''\ (4.35) 

For \p\>po, by (|42D, 

lko(t,p)|| + \\wi{t,p)\\ < Ce-^*e(p,0)||«;(0)||o. 
Using now the simple estimates 

e(j9,t)-i(e-^*p'/2^e-^*/2) < Ce(p,0)-\ 
for t > 1, e{p,t)-^ < Ce{p,Oy^ for t < 1, and 

for t > 1, m > 0, we conclude that 

||T,(t)-To(t)||i < C{l + t)-i\\w{0)\\o, (4.36) 
||^^Ki)-^oWlk < Cil + t)-'\\wiO)\\o. (4.37) 

The estimates (JSOlD and fl37[5|) follow now from Km . flCTj) . fICTD and fOTjl . □ 

Proof of Theorem 2. The proof goes along the lines of the one of Theorem 1 and we 
will be brief. We expand 

\c{W) = -\lw + -m(T, v) + n(T, v), 

where m{T,v) is the term quadradic in T but linear in v, and n{T,v) collects the terms 
that are quadratic and cubic in v. Note that, since LT = 0, -^Lw = ^Lv and thus 

DC{cu-^ + T)v = -Lv + m{T,v). (4.38) 

We write ( KT9^ as 

w(t) = exp( — -D^w(0) -\ — / dsexp( —D^)(m(s) + en(s))ds, 

e Jo e 



where 



D, = L + e^p ■ Vu, (4.39) 
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i.e. it equals (13.31) evaluated at ep. 

The decomposition w = T + ev yields, as in (13. 6p . (13. 7p : 

T{t) = Pe-'^^I''PT{Q) + tPe-'^^l''Qv{Q) + 
1 /"* 

- / c/sPe-(*-^)^^/^'g(m(s) + mis)) (4.40) 
e Jo 

v{t) = -Qe-*^^/^'PT(0) + Qe-*^^/^'Qt;(0) + 
e 

4 f dsQe-^'-'^''^/''Q(m(s) + en(s)) (4.41) 
Jo 

Consider first the case \p\ < po/e. We use Lemma 4.3, with t replaced by t/e^ and p 
by ep. This leads to 

T = T* + Ti (4.42) 

w = AT* + - / dsQQe-(*-')^'/^'QQm(s) + (4.43) 
e Jo 



where 



and 



T* = e-'P'^TiO) + [ dse-^'-'^P^^PBm{s), (4.44) 
Jo 



Ti = ePe-*^^/^'Q^7(0)+ f dsPe~^'-'^''^/'' Qn{s)+ f dsO{ep'e-''^'-''^P\\p\e-<'-''^''')m{s), 

Jo Jo 
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= Qe~'^^l' Qv{Q) + - / dsQe-^'-'^^^'' Qn{s) + e / dsO(|ppe-^(*-^)P )m(s). 

Then, following the proof of (iJS]), (lOTIl . ^^?M . (ICTl) . we get 

||ri(t,p)|| < Ce(l + t)--e(p,t)(l|r(0)|| + ||^;(0)||), (4.45) 
\\v,{t,p)\\ < C{t)e{l+t)-'e{p,mTm + ll^^(O)ll), (4.46) 

where C{t) = C max(t~^/^, log(l + 1)) < oo, for t > 0. We need to study the second term 
on the RHS of eq. (14. 43 p . Write m{s) = m{t) + (m(s) — m{t)) and consider the first term 

- / rfsQQe-(*-^)^^/^'QQm(t) = - / dsQQe-'''^/''QQm{t) + 0{e-'^/'') 
Jo Jo 

= QQD:'QQm{t) + 0(e-^*/^') 

= L-'m{t) + C(e) + 0(e-^*/^') (4.47) 

where we used QQD~^QQ = QL^^Q + (9(e) and where O(-) denotes a bound on the 
norm sup\p\^p^^/^e{p,t)^^\\-\\ of the remainders, since sup|p|<pg/^ e(p, )f:)~^||m(t)|| is bounded, 
uniformly in e, because of the bounds on T, v, coming from (14. 42p - (14.461) . 
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For the second term, we use the equations fl4.40p and (14.411) . to show that 
'^(i^)llt ^ C\t — s\, for \t — s\ small; whence 



\ f dsQQe-'^'-'^''^/''QQ{m{s) - m{t)) = 0{e). (4.48) 
e Jo 

Altogether we obtain 

v{t) = AT*{t) + L-^m{t) + V2 (4.49) 

and, for any t > 0, 

sup e(p,t)-^(||Ti|| + ||t;2||)^0 (4.50) 

|p|<Po/e 

as e — > 0. For |j9| > p^/e we have, from Proposition 4.2, 
which implies that 

sup e{p,tr\\\T\\ + \\v\\)^0, (4.51) 

|p|>poA 

as e ^ 0. Hence, in the space St, we have 

lim(T(t), vit)) = (T*(t), AT*{t) + L'^m{t)) = (T*(t), v*{t)) (4.52) 

Let us finally check that T* and v* satisfy the equations (K^ and (KT2\i . From fOSD 
and fl4.52p . we get 

AT* = -L-^DC{uo-^ + T*)v* (4.53) 

which is f l3.2ip if we recall the definition fl4.12p of A. 

From (14.441) and the definition (I4.12p of B, we see that T* satisfies 

t* = -p^kT* - —Pp ■ VujL^^m. (4.54) 
27r 

From (14.521) and the definition (I4.12p of A, we have 

L-^m = v* - AT* = v* + i{2n)-^L-^p ■ VuoT*. 

Substituting this into (I4.54p and recalling the definition of k in (13.111) . which implies that 
-p^nT* + p2(27r)-2p . VtuL-i . ^ equation (IX^ . □ 
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Proof of Lemma 4.3. We need to calculate Re~^^R' for R and R' either P or Q. We 
use the representation (14.141) for e~*^ and the formula (14.111) for P. We need the following 
expansions 

PP = P + PBQ + 0{p^), PP = P + AP + 0{p^) 

PQ = PBQ + Qp = Ap^ oip") 

PQ = P(l -P) = -PBQ + Oip^), QP = -QAP + 0{p^) 

QQ = {l-P)Q = Q- PBQ + QQ = Q_AP + 0{p^) 

where we used repeatedly PAP = {A maps E into E-^, because is odd in k). Thus 
we get, using the identities on the left, 

Pe~*^g = PPe-'^QQ + PQe'^^QQ 

= {P + PBQ)e-'^PBQ - PBQe-'^{Q - PBQ) + C(/e-^*f') 

= e-'P'^PBQ + 0{p^e~^'P') + 0{\p\e-''') (4.55) 

where we use the bounds of Proposition 4.2, B = 0{\p\), and Aq = fiap'^ + 0{p^), with 
lia being the eigenvalues of k, see (14. ip . 
Similarly, we get 



and 



Qe-*^Q = Ae-'P'^'PBQ + QQe-'^'QQ + 0{\pfe'''P') 



Pe-'^'P = e-'P"''P + Oi\p\e~''P) + 0{p^e-'') 
Qe-'^'P = Ae-'P'^P + Oip^e-^'^') + 0{\p\e-'''). 



□ 
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